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Abstract 

We discuss the theory of electromagnetic fields, with an emphasis on aspects 
relevant to radiofrequency systems in particle accelerators. We begin by re- 
viewing Maxwell's equations and their physical significance. We show that in 
free space, there are solutions to Maxwell's equations representing the propa- 
gation of electromagnetic fields as waves. We introduce electromagnetic po- 
tentials, and show how they can be used to simplify the calculation of the fields 
in the presence of sources. We derive Poynting's theorem, which leads to ex- 
pressions for the energy density and energy flux in an electromagnetic field. 
We discuss the properties of electromagnetic waves in cavities, waveguides 
and transmission lines. 

1 Maxwell's equations 

Maxwell's equations may be written in differential form as follows: 

(1) 
(2) 

(3) 
(4) 

determine the force on a particle of 



The electric displacement D and magnetic intensity H are related to the electric field and magnetic flux 
density by the constitutive relations: 

D = eE, 
B = nH. 

The electric permittivity e and magnetic permeability depend on the medium within which the fields 
exist. The values of these quantities in vacuum are fundamental physical constants. In SI units: 

/xo = 47r X 10"'^Hm-\ 
1 



where c is the speed of light in vacuum. The permittivity and permeability of a material characterize the 
response of that material to electric and magnetic fields. In simplified models, they are often regarded 
as constants for a given material; however, in reality the permittivity and permeability can have a com- 
plicated dependence on the fields that are present. Note that the relative permittivity and the relative 
permeability fir are frequently used. These are dimensionless quantities, defined by: 



V-D = p, 
VB = 0, 

„ DB 

The fields B (magnetic flux density) and E (electric field strength) 
charge q travelling with velocity v (the Lorentz force equation): 



F = q(E + vy.B \ . 
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Fig. 1: Snapshot of a numerical solution to Maxwell's equations for a bunch of electrons moving through a beam 
position monitor in an accelerator vacuum chamber The colours show the strength of the electric field. The bunch 
is moving from right to left: the location of the bunch corresponds to the large region of high field intensity towards 
the left hand side. (Image courtesy of M. Korostelev.) 

That is, the relative permittivity is the permittivity of a material relative to the permittivity of free space, 
and similarly for the relative permeability. 

The quantities p and J are respectively the electric charge density (charge per unit volume) and 
electric current density ( J • n is the charge crossing unit area perpendicular to unit vector ft per unit time). 
Equations ^ and are independent of p and J, and are generally referred to as the "homogeneous" 
equations; the other two equations, ([l]) and (jsjl are dependent on p and J, and are generally referred to 
as the "inhomogeneous" equations. The charge density and current density may be regarded as sources 
of electromagnetic fields. When the charge density and current density are specified (as functions of 
space, and, generally, time), one can integrate Maxwell's equations ([l]l-(|3]) to find possible electric and 
magnetic fields in the system. Usually, however, the solution one finds by integration is not unique: for 
example, as we shall see, there are many possible field patterns that may exist in a cavity (or waveguide) 
of given geometry. 

Most realistic situations are sufficiently complicated that solutions to Maxwell's equations cannot 
be obtained analytically. A variety of computer codes exist to provide solutions numerically, once the 
charges, currents, and properties of the materials present are all specified, see for example References 
Solving for the fields in realistic systems (with three spatial dimensions, and a dependence on 
time) often requires a considerable amount of computing power; some sophisticated techniques have 
been developed for solving Maxwell's equations numerically with good efficiency Q. An example of a 
numerical solution to Maxwell's equations in the context of a particle accelerator is shown in Fig.[T] We 
do not consider such techniques here, but focus instead on the analytical solutions that may be obtained 
in idealized situations. Although the solutions in such cases may not be sufficiently accurate to complete 
the design of real accelerator components, the analytical solutions do provide a useful basis for describing 
the fields in (for example) real RF cavities and waveguides. 

An important feature of Maxwell's equations is that, for systems containing materials with con- 
stant permittivity and permeability (i.e. permittivity and permeability that aie independent of the fields 
present), the equations are linear in the fields and sources. That is, each term in the equations involves 
a field or a source to (at most) the first power, and products of fields or sources do not appear. As a 
consequence, the principle of superposition applies: if Ei,Bi and E2,B2 are solutions of Maxwell's 
equations with given boundary conditions, then = Ei + E2 and = Bi + B2 will also be so- 
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lutions of Maxwell's equations, with the same boundary conditions. This means that it is possible to 
represent complicated fields as superpositions of simpler fields. An important and widely used analysis 
technique for electromagnetic systems, including RF cavities and waveguides, is to find a set of solu- 
tions to Maxwell's equations from which more complete and complicated solutions may be constructed. 
The members of the set are known as modes; the modes can generally be labelled using mode indices. 
For example, plane electromagnetic waves in free space may be labelled using the three components of 
the wave vector that describes the direction and wavelength of the wave. Important properties of the 
electromagnetic fields, such as the frequency of oscillation, can often be expressed in terms of the mode 
indices. 

Solutions to Maxwell's equations lead to a rich diversity of phenomena, including the fields around 
charges and currents in certain basic configurations, and the generation, transmission and absorption of 
electromagnetic radiation. Many existing texts cover these phenomena in detail; for example. Grant 
and Phillips [5], or the authoritative text by Jackson [6]. We consider these aspects rather briefly, with 
an emphasis on those features of the theory that are important for understanding the properties of RF 
components in accelerators. 

2 Integral theorems and the physical interpretation of Maxwell's equations 
2.1 Gauss' theorem and Coulomb's law 

Guass' theorem states that for any smooth vector field a: 

[ V -adV = I a-dS, 

Jv JdV 

where F is a volume bounded by the closed surface dV. Note that the area element dS is oriented to 
point out of V. 

Gauss' theorem is helpful for obtaining physical interpretations of two of Maxwell's equations, 
([T]l and ([2]). First, applying Gauss' theorem to ([T]l gives: 

[ V -DdV = (f D-dS = q, (6) 

Jv JdV 

where q = Jy p dV is the total charge enclosed by dV. 

Suppose that we have a single isolated point charge in an homogeneous, isotropic medium with 
constant permittivity e. In this case, it is interesting to take dV to be a sphere of radius r. By symmetry, 
the magnitude of the electric field must be the same at all points on dV, and must be normal to the 
surface at each point. Then, we can perform the surface integral in ([6]): 

/ D ■ dS = Anr'^D. 
JdV 

This is illustrated in Fig.|2} the outer circle represents a cross-section of a sphere (dV) enclosing volume 
V, with the charge q at its centre. The red arrows in Fig. |2|represent the electric field lines, which are 
everywhere perpendicular to the surface dV. Since D = eE, we find Coulomb's law for the magnitude 
of the electric field around a point charge: 

E = ^—. 

Applied to Maxwell's equation Q, Gauss' theorem leads to: 

/ V -BdV = I B-dS = 0. 

Jv JdV 
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In other words, the magnetic flux integrated over any closed surface must equal zero - at least, until we 
discover magnetic monopoles. Lines of magnetic flux always occur in closed loops; lines of electric field 
may occur in closed loops, but in the presence of electric charges will have start (and end) points on the 
electric charges. 

2.2 Stokes' theorem, Ampere's law, and Faraday's law 

Stokes' theorem states that for any smooth vector field a: 



V xa-dS = f a-dl, (7) 
S JdS 

where the closed loop dS bounds the surface S. Applied to Maxwell's equation ([3]), Stokes' theorem 
leads to: 



H-dl= / J-dS, (8) 
as Js 

which is Ampere's law. From Ampere's law, we can derive an expression for the strength of the magnetic 
field around a long, straight wire carrying current /. The magnetic field must have rotational symmetry 
around the wire. There are two possibilities: a radial field, or a field consisting of closed concentric 
loops centred on the wire (or some superposition of these fields). A radial field would violate Maxwell's 
equation Q. Therefore, the field must consist of closed concentric loops; and by considering a circular 
loop of radius r, we can perform the integral in Eq. ([8]): 

2TTrH = I, 

where / is the total current carried by the wire. In this case, the line integral is performed around a loop 
dS centered on the wire, and in a plane perpendicular to the wire: essentially, this corresponds to one of 
the magnetic field lines, see Fig. [3] The total current passing through the surface S bounded by the loop 
dS is simply the total current /. 

In an homogeneous, isotropic medium with constant permeability fi, B = fioH, and we obtain the 
expression for the magnetic flux density at distance r from the wire: 

B = i^- (9) 
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Fig. 3: Magnetic field lines around a long straight wire carrying a current I. 



Finally, applying Stokes' theorem to the homogeneous Maxwell's equation Q, we find: 

E-dl = -— I B-dS. (10) 



Ids 9t js 

Defining the electromotive force S' as the integral of the electric field around a closed loop, and the 



magnetic flux $ as the integral of the magnetic flux density over the surface bounded by the loop, Eq. ( 10 1 
gives: 

which is Faraday's law of electromagnetic induction. 

Maxwell's equations Q and Q are significant for RF systems: they tell us that a time dependent 
electric field will induce a magnetic field; and a time dependent magnetic field will induce an electric 
field. Consequently, the fields in RF cavities and waveguides always consist of both electric and magnetic 
fields. 



3 Electromagnetic waves in free space 

In free space (i.e. in the absence of any charges or currents) Maxwell's equations have a trivial solution 
in which all the fields vanish. However, there are also non-trivial solutions with considerable practical 
importance. In general, it is difficult to write down solutions to Maxwell's equations, because two of the 
equations involve both the electric and magnetic fields. However, by taking additional derivatives, it is 
possible to write equations for the fields that involve only either the electric or the magnetic field. This 
makes it easier to write down solutions: however, the drawback is that instead of first-order differential 
equations, the new equations are second-order in the derivatives. There is no guarantee that a solution 
to the second-order equations will also satisfy the first-order equations, and it is necessary to impose 
additional constraints to ensure that the first-order equations are satisfied. Fortunately, it turns out that 
this is not difficult to do, and taking additional derivatives is a useful technique for simplifying the 
analytical solution of Maxwell's equations in simple cases. 



3.1 Wave equation for the electric field 

In free space. Maxwell's equations - Q take the form: 



V-E = 0, (12) 
V-B = 0, (13) 
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V X B 

V X E 



1 dE 
dB 



where we have defined: 



(14) 
(15) 

(16) 



Our goal is to find a form of the equations in which the fields E and B appear separately, and not together 



in the same equation. As a first step, we take the curl of both sides of Eq. ( 15 1, and interchange the order 
of differentiation on the right hand side (which we are allowed to do, since the space and time coordinates 
are independent). We obtain: 

d 

\7 xV X E = -—V X B. (17) 
ot 



Substituting for V x i? from Eq. ( 14 1, this becomes: 



V X V X £; 



1 d^E 



(18) 



This second-order differential equation involves only the electric field, E, so we have achieved our 
aim of decoupling the field equations. However, it is possible to make a further simplification, using a 
mathematical identity. For any differentiable vector field a: 



V xV X a = V(V • a) - V^a. 



Using the identity ( [19] ), and also making use of Eq. ( [12] ), we obtain finally: 

1 d'^E 



V'^E 



c2 df^ 



0. 



(19) 



(20) 



Eq. ( [20] ) is the wave equation in three spatial dimensions. Note that each component of the electric field 
independently satisfies the wave equation. The solution, representing a plane wave propagating in the 
direction of the vector k, may be written in the form: 



E = Eq cos (k ■ f — Ult + 



(21) 



where: Eq is a constant vector; (/}q is a constant phase; uj and k are constants related to the frequency / 
and wavelength A of the wave by: 



u = 27r/, 
27r 



A 



1^1 



(22) 
(23) 



If we substitute Eq. ( [21) into the wave equation ( \20} , we find that it provides a valid solution as long as 
the angular frequency uj and wave vector k satisfy the dispersion relation: 



\k\ 



c. 



(24) 



If we inspect Eq. (21 ), we see that a particle travelling in the direction of k has to move at a speed 
io/\k\ in order to remain at the same phase in the wave: thus, the quantity c is the phase velocity of 
the wave. This quantity c is, of course, the speed of light in a vacuum; and the identification of light 
with an electromagnetic wave (with the phase velocity related to the electric permittivity and magnetic 
permeabihty by Eq. ( 16 1) was one of the great achievements of 19"^ century physics. 
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3.2 Wave equation for the magnetic field 

So far, we have only considered the electric field. But Maxwell's equation Q tells us that an electric 
field that varies with time must have a magnetic field associated with it. Therefore, we should look for a 



(non-trivial) solution for the magnetic field in free space. Starting with Eq. ( 14 1, and following the same 
procedure as above, we find that the magnetic field also satisfies the wave equation: 



1 d'^B 



0, 



with a similar solution: 



B = Bq cos 



i(k- 



ijjt + 



(25) 
(26) 



Here, we have written the same constants w, k and 0o as we used for the electric field, though we do 
not so far know they have to be the same. We shall show in the following section that these constants do 
indeed need to be the same for both the electric field and the magnetic field. 



3.3 Relations between electric and magnetic fields in a plane wave in free space 

As we commented above, although taking additional derivatives of Maxwell's equations allows us to 
decouple the equations for the electric and magnetic fields, we must impose additional constraints on the 
solutions, to ensure that the first-order equations are satisfied. In particular, substituting the expressions 



for the fields (21 1 and (26 1 into Eqs. ( 12 1 and ( 13 1 respectively, and noting that the latter equations must 
be satisfied at all points in space and at all times, we obtain: 



k-EQ 
k-Bo 



0, 
0. 



(27) 
(28) 



Since k represents the direction of propagation of the wave, we see that the electric and magnetic fields 
must at all times and all places be perpendicular to the direction in which the wave is travelling. This is 
a feature that does not appear if we only consider the second-order equations. 



Finally, substituting the expressions for the fields (21 1 and (26) into Eqs. ( 15 I and ( 14 1 respectively. 



and again noting that the latter equations must be satisfied at all points in space and at all times, we see 
first that the quantities uj, k and appearing in (21 1 and (26 1 must be the same in each case. Also, we 
have the following relations between the magnitudes and directions of the fields: 



kyi E< 







ujBn 



(29) 
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kx Bo = -ojEq. (30) 

If we choose a coordinate system so that Eq is parallel to the x axis and Bq is parallel to the y axis, 
then k must be parallel to the z axis: note that the vector product ^ x i3 is in the same direction as the 
direction of propagation of the wave - see Fig.|4] The magnitudes of the electric and magnetic fields are 
related by: 

Note that the wave vector k can be chosen arbitrarily: there are infinitely many "modes" in which 
an electromagnetic wave propagating in free space may appear; and the most general solution will be 
a sum over all modes. When the mode is specified (by giving the components of k), the frequency 



is determined from the dispersion relation (24i. However, the amplitude and phase are not determined 



(although the electric and magnetic fields must have the same phase, and their amplitudes must be related 



by Eq. 

Finally, note that all the results derived in this section are strictly true only for electromagnetic 
fields in a vacuum. The generalisation to fields in uniform, homogenous, linear (i.e. constant perme- 
ability /i and permittivity e) nonconducting media is straightforward. However, new features appear for 
waves in conductors, on boundaries, or in nonlinear media. 



3.4 Complex notation for electromagnetic waves 

We finish this section by introducing the complex notation for free waves. Note that the electric field 
given by equation ( pT) can also be written as: 

^ = Re ^oe*'^°e^(^-^^-'^*\ (32) 

To avoid continually writing a constant phase factor when dealing with solutions to the wave equation, 
we replace the real (constant) vector Eq by the complex (constant) vector Eq = Eqc^'^^. Also, we note 
that since all the equations describing the fields are linear, and that any two solutions can be linearly 
superposed to construct a third solution, the complex vectors: 

E' = ^e*^^'^^"''*^ (33) 
B' = S^e*(^-^^-^*) (34) 

provide mathematically valid solutions to Maxwell's equations in free space, with the same relationships 
between the various quantities (frequency, wave vector, amplitudes, phase) as the solutions given in 



Eqs. (21 1 and (26 1. Therefore, as long as we deal with linear equations, we can carry out all the algebraic 
manipulation using complex field vectors, where it is implicit that the physical quantities are obtained 
by taking the real parts of the complex vectors. However, when using the complex notation, particular 
care is needed when taking the product of two complex vectors: to be safe, one should always take the 
real part before multiplying two complex quantities, the real parts of which represent physical quantities. 
Products of the electromagnetic field vectors occur in expressions for the energy density and energy flux 
in an electromagnetic field, as we shall see below. 



4 Electromagnetic waves in conductors 

Electromagnetic waves in free space are characterized by an amplitude that remains constant in space and 
time. This is also true for waves travelling through any isotropic, homogeneous, linear, non-conducting 
medium, which we may refer to as an "ideal" dielectric. The fact that real materials contain electric 
charges that can respond to electromagnetic fields means that the vacuum is really the only ideal dielec- 
tric. Some real materials (for example, many gases, and materials such as glass) have properties that 
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approximate those of an ideal dielectric, at least over certain frequency ranges: such materials are trans- 
parent. However, we know that many materials are not transparent: even a thin sheet of a good conductor 
such as aluminium or copper, for example, can provide an effective barrier for electromagnetic radiation 
over a wide range of frequencies. 

To understand the shielding effect of good conductors is relatively straightforward. Essentially, we 
follow the same procedure to derive the wave equations for the electromagnetic fields as we did for the 
case of a vacuum, but we include additional terms to represent the conductivity of the medium. These 
additional terms have the consequence that the amplitude of the wave decays as the wave propagates 
through the medium. The rate of decay of the wave is characterised by the skin depth, which depends 
(amongst other things) on the conductivity of the medium. 

Let us consider an ohmic conductor. An ohmic conductor is defined by the relationship between 
the current density J at a point in the conductor, and the electric field E existing at the same point in the 
conductor: 

J = aE, (35) 

where o" is a constant, the conductivity of the material. 

In an uncharged ohmic conductor. Maxwell's equations Q - Q take the form: 

V-E = 0, (36) 
V-B = 0, (37) 
dE 

VxB = uaE + Lie—, (38) 

ot 

dB 

VxE = --, ,39) 

where ^ is the (absolute) permeability of the medium, and e is the (absolute) permittivity. Notice the 
appearance of the additional term on the right hand side of Eq. p8|), compared to Eq. ([14]). 



Following the same procedure as led to Eq. (20 1, we derive the following equation for the electric 
field in a conducting medium: 

^ - dE d'^E 
V^E-,a--,e^ = 0. (40) 

This is again a wave equation, but with a term that includes a first-order time derivative. In the equation 
for a simple harmonic oscillator, such a term would represent a "frictional" force that leads to dissipation 
of the energy in the oscillator. There is a similar effect here; to see this, let us try a solution of the same 
form as for a wave in free space. The results we are seeking can be obtained more directly if we used the 
complex notation: 

^ = ^oe*(^-^-^*). (41) 



Substituting into the wave equation (40 1, we obtain the dispersion relation: 



-P + i^auj + fj.euj'^ = 0. (42) 



Let us assume that the frequency w is a real number. Then, to find a solution to Eq. ( [42| l, we have to 
allow the wave vector k to be complex. Let us write the real and imaginary parts as a and (3 respectively: 

k = a + i(3. (43) 



Substituting (43 1 into (42 1 and equating real and imaginary parts, we find (after some algebra) that: 
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X 




Fig. 5: Electric and magnetic fields in a plane electromagnetic wave in a conductor. The wave vector is in the 
direction of the +z axis. 



I?l = f5[. (45) 

2\a\ 

To understand the physical significance of a. and /3, we write the solution (41 1 to the wave equation as: 

^ = ^oe-'^-V("-^^-^*). (46) 

We see that there is still a wave-like oscillation of the electric field, but there is now also an exponential 
decay of the amplitude. The wavelength is determined by the real part of the wave vector: 

X=%. (47) 

|q| 

The imaginary part of the wave vector gives the distance 5 over which the amplitude of the wave falls by 
a factor 1/e, known as the skin depth: 

(5 = i. (48) 

1/31 

Accompanying the electric field, there must be a magnetic field: 

B = Boe'^'^-'^-'^^l (49) 

From Maxwell's equation Q, the amplitudes of the electric and magnetic fields must be related by: 

kx Eo = ujBq. (50) 

The electric and magnetic fields are perpendicular to each other, and to the wave vector: this is the 
same situation as occurred for a plane wave in free space. However, since k is complex for a wave in 
a conductor, there is a phase difference between the electric and magnetic fields, given by the complex 
phase of k. The fields in a plane wave in a conductor are illustrated in Fig. pi 



The dispersion relation ( 42 1 gives a rather complicated algebraic relationship between the fre- 
quency and the wave vector, in which the electromagnetic properties of the medium (permittivity, per- 
meabihty and conductivity) all appear. However, in many cases it is possible to write much simpler 
expressions that provide good approximations. First, there is the "poor conductor" regime: 

if cr ^ we, then I a| a; ^/Tte, |/3| ~ ^W— . (51) 
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The wavelength is related to the frequency in the way that we would expect for a dielectric. 
Next there is the "good conductor" regime: 



if cr » ue, then |a| ~ y ~2~' '^^'^ ^ 

Here the situation is very different. The wavelength depends directly on the conductivity: for a good 
conductor, the wavelength is very much shorter than it would be for a wave at the same frequency in free 
space. The real and imaginary parts of the wave vector are approximately equal: this means that there 
is a significant reduction in the amplitude of the wave even over one wavelength. Also, the electric and 
magnetic fields are approximately tt/A out of phase. 

The reduction in amplitude of a wave as it travels through a conductor is not difficult to understand. 
The electric charges in the conductor move in response to the electric field in the wave. The motion of 
the charges constitutes an electric current in the conductor, which results in ohmic losses: ultimately, the 
energy in the wave is dissipated as heat in the conductor. Note that whether or not a given material can be 
described as a "good conductor" depends on the frequency of the wave (and permittivity of the material): 
at a high enough frequency, any material will become a poor conductor. 



5 Energy in electromagnetic fields 

Waves are generally associated with the propagation of energy: the question then arises as to whether 
this is the case with electromagnetic waves, and, if so, how much energy is carried by a wave of a given 
amplitude. To address this question, we first need to find general expressions for the energy density and 
energy flux in an electromagnetic field. The appropriate expressions follow from Poynting's theorem, 
which may be derived from Maxwell's equations. 



5.1 Poynting's theorem 

First, we take the scalar product of Maxwell's equation Q with the magnetic intensity H on both sides, 
to give: 

- dB 

H -V X E = -H ■ —. (53) 
ot 

Then, we take the scalar product of (jsjl with the electric field E on both sides to give: 

E-V X H = E- J + E-^. (54) 

ot 

Now we subtract Eq. (|54]) from Eq. ([53]) to give: 



H-VxE-E-VxH = -E-J-E- — -H-—. (55) 

dt dt 



This may be rewritten as: 



^J^eE^ + ^liH'] =-V-(ExH)-E-J. (56) 



Equation (56 1 is Poynting's theorem. It does not appear immediately to tell us much about the energy 
in an electromagnetic field; but the physical interpretation becomes a little clearer if we convert it from 
differential form into integral form. Integrating each term on either side over a volume V, and changing 
the first term on the right hand side into an integral over the closed surface A bounding V, we write: 



dt 



[ (Ue + Uh) dV = - (f S-dA- [ E-JdV, (57) 

Jv J A Jv 
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where: 



Ue = \eE'' (58) 

Uh = \liH^ (59) 
S = Ex H. (60) 



The physical interpretation follows from the volume integral on the right hand side of Eq. (57i: 
this represents the rate at which the electric field does work on the charges contained within the volume 
V. If the field does work on the charges within the field, then there must be energy contained within 
the field that decreases as a result of the field doing work. Each of the terms within the integral on the 



left hand side of Eq. (57 1 has the dimensions of energy density (energy per unit volume). Therefore, the 
integral has the dimensions of energy; it is then natural to interpret the full expression on the left hand 
side of Eq. ( [57| ) as the rate of change of energy in the electromagnetic field within the volume V. The 
quantities U e and Uh represent the energy per unit volume in the electric field and in the magnetic field 
respectively. 



Finally, there remains the interpretation of the first term on the right hand side of Eq. (57i. As 
well as the energy in the field decreasing as a result of the field doing work on charges, the energy may 
change as a result of a flow of energy purely within the field itself (i.e. even in the absence of any 



electric charge). Since the first term on the right hand side of Eq. (57 1 is a surface integral, it is natural 



to interpret the vector inside the integral as the energy flux within the field, i.e. the energy crossing unit 



area (perpendicular to the vector) per unit time. The vector S defined by Eq. (60 1 is called the Poynting 
vector. 



5.2 Energy in an electromagnetic wave 

As an application of Poynting's theorem (or rather, of the expressions for energy density and energy flux 
that arise from it), let us consider the energy in a plane electromagnetic wave in free space. As we noted 
above, if we use complex notation for the fields, then we should take the real part to find the physical 
fields before using expressions involving the products of fields (such as the expressions for the energy 
density and energy flux). 

The electric field in a plane wave in free space is given by: 



E = Eocos(k-f-ujt + (j)o)- (61) 
Thus, the energy density in the electric field is: 



Ue = ^SqE"^ = -eqEI cos^ [k ■ r - ujt + 4>q] . (62) 



1 1 

-eo^ = -t 

2 2 

If we take the average over time at any point in space (or, the average over space at any point in time), 
we find that the average energy density is: 

{Ue) = \eoEl (63) 
The magnetic field in a plane wave in free space is given by: 

B = Bocos[k ■ r - ojt + 0oj , (64) 

where: 

l^ol = ^. (65) 
c 
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Since B = hqH, the energy density in the magnetic field is: 

Uh = ^i^qH^ "2^^° (k ■ f - ujt + (l)oy (66) 

If we take the average over time at any point in space (or, the average over space at any point in time), 
we find that the average energy density is: 

{Uh) = (67) 
4/io 

Using the relationship ([65|) between the electric and magnetic fields in a plane wave, this can be written: 



iUH) = ^^- (68) 



Then, using Eq. ( 16 1: 



{Uh) = leoEl (69) 



We see that in a plane electromagnetic wave in free space, the energy is shared equally between 
the electric field and the magnetic field, with the energy density averaged over time (or, over space) given 
by: 

{U) = ^eoEl (70) 



Finally, let us calculate the energy flux in the wave. For this, we use the Poynting vector (6O1: 

S = ExH = k—E^cos'^(k-r-ujt + 4>o), (71) 

where is a unit vector in the direction of the wave vector. The average value (over time at a particular 
point in space, or over space at a particular time) is then given by: 

{S) = l—Efk = leocEfk = ^k, (72) 

where Zq is the impedance of free space, defined by: 

Zo = ^[^. (73) 



Zq is a physical constant, with value Zq w 376.73 U. Using Eq. (70 1 we find the relation between energy 
flux and energy density in a plane electromagnetic wave in free space: 

{S) = {U)ck. (74) 

This is the relationship that we might expect in this case: the mean energy flux is given simply by the 
mean energy density moving at the speed of the wave in the direction of the wave. But note that this is 
not the general case. More generally, the energy in a wave propagates at the group velocity, which may 
be different from the phase velocity. For a plane electromagnetic wave in free space, the group velocity 
happens to be equal to the phase velocity, c. We shall discuss this further when we consider energy 
propagation in waveguides. 
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6 Electromagnetic potentials 

We have seen that to find non-trivial solutions for Maxwell's electromagnetic field equations in free 
space, it is helpful to take additional derivatives of the equations, since this allows us to construct separate 
equations for the electric and magnetic fields. The same technique can be used to find solutions for the 
fields when sources (charge densities and currents) are present. Such situations are important, since they 
arise in the generation of electromagnetic waves. However, it turns out that in systems where charges 
and currents are present, it is often simpler to work with the electromagnetic potentials, from which the 
fields may be obtained by differentiation, than with the fields directly. 

6.1 Relationships between the potentials and the fields 

The scalar potential (j) and vector potential A are defined so that the electric and magnetic fields are 
obtained using the relations: 

B = V X 1, (75) 
BA 

E = -V0-^. (76) 
We shall show below that as long as (f) and A satisfy appropriate equations, then the fields B and E 



derived from them using Eqs. ( |75| ) and ( |76| ) satisfy Maxwell's equations. But first, note that there is a 
many-to-one relationship between the potentials and the fields. That is, there are many different poten- 
tials that can give the same fields. For example, we could add any uniform (independent of position) 
value to the scalar potential 0, and leave the electric field E unchanged, since the gradient of a constant 
is zero. Similarly, we could add any vector function with vanishing curl to the vector potential A; and 
if this function is independent of time, then again the electric and magnetic fields are unchanged. This 
property of the potentials is known as gauge invariance, and is of considerable practical value, as we 
shall see below. 



6.2 Equations for the potentials 

The fact that there is a relationship between the potentials and the fields implies that the potentials that 
are allowed in physics have to satisfy certain equations, corresponding to Maxwell's equations. This is, 
of course, the case. In this section, we shall derive the equations that must be satisfied by the potentials, 
if the fields that are derived from them are to satisfy Maxwell's equations. 

However, to begin with, we show that two of Maxwell's equations (the ones independent of the 



sources) are in fact satisfied if the fields are derived from any potentials (p and A using Eqs. (75) and 



(76 1. First, since the divergence of the curl of any differentiable vector field is always zero: 

V • V X i* = 0, (77) 

it follows that Maxwell's equation ^ is satisfied for any vector potential A. Then, since the curl of the 
gradient of any differentiable scalar field is always zero: 

V X V(/> = 0, (78) 

Maxwell's equation (Q is satisfied for any scalar potential (p and vector potential A (as long as the 
magnetic field is obtained from the vector potential by Eq. ([75])). 

Now let us consider the equations involving the source terms (the charge density p and current 



density J). Differential equations for the potentials can be obtained by substituting from Eqs. (75 1 
and ( [76] ) into Maxwell's equations ([T]) and ([3]). We also need to use the constitutive relations between 
the magnetic field B and the magnetic intensity H, and between the electric field E and the electric 
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displacement D. For simplicity, let us assume a system of charges and currents in free space; then the 
constitutive relations are: 

D = eoE, B = fioH. (79) 



Substituting from Eq. (76 1 into Maxwell's equation ([T]) gives: 



ot eo 



Similarly, substituting from Eq. ( [75| ) into Maxwell's equation Q gives (after some algebra): 



(80) 



(81) 



Eqs. (80 1 and (81 1 relate the electromagnetic potentials to a charge density p and current density J in 
free space. Unfortunately, they are coupled equations (the scalar potential (p and vector potental J each 
appear in both equations), and are rather complicated. However, we noted above that the potentials for 
given electric and magnetic fields are not unique: the potentials have the property of gauge invariance. By 
imposing an additional constraint on the potentials, known as a gauge condition, it is possible to restrict 
the choice of potentials. With an appropriate choice of gauge, it is possible to decouple the potentials, 
and furthermore, arrive at equations that have standard solutions. In particular, with the gauge condition: 



V-A + 



1 



dt 



then Eq. ( 80 1 becomes: 



and Eq. ( 81 1 becomes: 



1 d^ct) 



1 d^A 



0, 

_P_ 



(82) 



(83) 



(84) 



Eqs. ( [83] ) and ( 84 1 have the form of wave equations with source terms. It is possible to write solutions in 
terms of integrals over the sources: we shall do this shortly. However, before we do so, it is important to 
note that for any given potentials, it is possible to find new potentials that satisfy Eq. ([82]), but give the 



same fields as the original potentials. Eq. ( 82 1 is called the Lorenz gauge. The proof proceeds as follows. 

First we show that any function ijj of position and time can be used to construct a gauge trans- 
formation; that is, we can use ijj to find new scalar and vector potentials (different from the original 
potentials) that given the same electric and magnetic fields as the original potentials. Given the original 
potentials (j) and A, and a function ip, let us define new potentials 0' and A' , as follows: 



+ 



A' 



A 



'dt' 



(85) 
(86) 



Eqs. ( 85 1 and ( 86 1 represent a gauge transformation. If the original potentials give fields E and B, then 
the magnetic field derived from the new vector potential is: 



B' = \7 X A' = V X A = B, 



(87) 



where we have used the fact that the curl of the gradient of any scalar function is zero. The electric field 
derived from the new potentials is: 



E' 



dA' 
~dt 



15 



(88) 



Here, we have made use of the fact that position and time are independent variables, so it is possible to 
interchange the order of differentiation. We see that for any function if:, the fields derived from the new 
potentials are the same as the fields derived from the original potentials. We say that ip generates a gauge 
transformation: it gives us new potentials, while leaving the fields unchanged. 

Finally, we show how to choose a gauge transformation so that the new potentials satisfy the 
Lorenz gauge condition. In general, the new potentials satisfy the equation: 

Suppose we have potentials (j) and A that satsify: 

where / is some function of position and time. (If / is non-zero, then the potentials 4> and A do not 
satisfy the Lorenz gauge condition.) Therefore, if ■0 satisfies: 

Then the new potentials (/)' and A' satisfy the Lorenz gauge condition: 

c"' ot 



Notice that Eq. (91 1 again has the form of a wave equation, with a source term. Assuming that we can 
solve such an equation, then it is always possible to find a gauge transformation such that, starting from 
some given original potentials, the new potentials satisfy the Lorenz gauge condition. 

6.3 Solution of the wave equation with source term 

^ r del) 



In the Lorenz gauge ( 82 1 



the vector potential A and the scalar potential cp satisfy the wave equations (84i and (83 1 : 



15^0 _ P 



Note that the wave equations have the form (for given charge density and current density) of two uncou- 
pled second-order differential equations. In many situations, it is easier to solve these equations, than 
to solve Maxwell's equations for the fields (which take the form of four first-order coupled differential 
equations). 
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origin 

Fig. 6: Integration over a distributed source. 



We state (without giving a full proof) the solution to the wave equation with a source term. The 
solution is: 

where: 

\r' — r\ 



t' = t 



(94) 



The integral extends over all space - see Fig. [6] Note that the source at each point in the integral has to be 
evaluated at a time t', which depends on the distance between the source point and the obeservation point 
(at which we are evaluating the potential). From the equation in the absence of any sources, we expect 
variations in the potential to propagate through space at a speed c (the speed of light). The difference 
between t' and t simply accounts for the time taken for the effect of any change in the charge density at 
the source point to propagate through space to the observation point. 



We do not present a proof that Eq. ([93) represents a solution to the wave equation (83 1. However, 



we can at least see that in the static case, Eq. (83 1 reduces to the familiar form of Poisson's equation: 



V 



2^ ^ _P_ 



For a point charge q ata point ro, Eq. (95 1 has solution: 

1 q 



Atteo \ro-f\' 

This can be obtained directly from Eq. ( [93] ) if the source is given by a Dirac delta function: 

p(f ) = qS^f' - fo). 



(95) 



(96) 



(97) 



The wave equation for the vector potential 
form to that for the scalar potential: 



has a solution that can be expressed in a similar 



(98) 



Note that Eq. ([98]) can be expressed as three independent equations (integrals), for the three components 
of the vector potential. 
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6.4 Physical significance of the fields and potentials 

An electromagnetic field is really a way of describing tiie interaction between particles that have elec- 
tric charges. Given a system of charged particles, one could, in principle, write down equations for the 
evolution of the system purely in terms of the positions, velocities, and charges of the various particles. 
However, it is often convenient to carry out an intermediate step in which one computes the fields gener- 
ated by the particles, and then computes the effects of the fields on the motion of the particles. Maxwell's 
equations provide a prescription for computing the fields arising from a given system of charges. The 
effects of the fields on a charged particle are expressed by the Lorentz force equation: 

E + vxbY (99) 

where F is the force on the particle, q is the charge, and v is the velocity of the particle. The motion of 
the particle under the influence of a force F is then given by Newton's second law of motion: 

—-ymv = F. (100) 
at 



Eqs. (99 1 and (100 1 make clear the physical significance of the fields. But what is the significance of 
the potentials? At first, the feature of gauge invariance appears to make it difficult to assign any definite 
physical significance to the potentials: in any given system, we have a certain amount of freedom in 
changing the potentials without changing the fields that are present. However, let us consider first the 
case of a particle in a static electric field. In this case, the Lorentz force is given by: 

F = qE = -qV(l). (101) 

If the particle moves from position fi to position 1^2 under the influence of the Lorentz force, then the 
work done on the particle (by the field) is: 

W = f^F- di= -q / ' V0 • de = -q [(/.(fa) - ^{n)] . (102) 

J ri J ri 

Note that the work done by the field when the particle moves between two points depends on the dif- 
ference in the potential at the two points; and that the work done is independent of the path taken by 
the particle in moving between the two points. This suggests that the scalar potential (p is related to the 
energy of a particle in an electrostatic field. If a time-dependent magnetic field is present, the analysis 
becomes more complicated. 

A more complete understanding of the physical significance of the scalar and vector potentials 
is probably best obtained in the context of Hamiltonian mechanics. In this formalism, the equations of 
motion of a particle are obtained from the Hamiltonian, H{x,p; t); the Hamiltonian is a function of the 
particle coordinates x, the (canonical) momentum p, and an independent variable t (often corresponding 
to the time). Note that the canonical momentum can (and generally does) differ from the usual mechani- 
cal momentum. The Hamiltonian defines the dynamics of a system, in the same way that a force defines 
the dynamics in Newtonian mechanics. In Hamiltonian mechanics, the equations of motion of a particle 
are given by Hamilton's equations: 

dxi dH (103) 

dt dpi ' 

dpi dH (104) 

dt dxi 

In the case of a charged particle in an electromagnetic field, the Hamiltonian is given by: 



H = cJ{p-qA)^+m^c'^+q(j), (105) 
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(106) 



where the canonical momentum is: 

p = (3^mc + qA, 

where /3 is the normalised velocity of the particle, P = v/c. 

Note that Eqs. ( |103| ) - ( |106| ) give the same dynamics as the Lorentz force equation ( 99 1 together 



with Newton's second law of motion, Eq. (100 1; they are just written in a different formalism. The 
significant point is that in the Hamiltonian formalism, the dynamics are expressed in terms of the po- 
tentials, rather than the fields. The Hamiltonian can be interpreted as the "total energy" of a particle, £. 



Combining equations ( 105 1 and ( 106 1, we find: 



(107) 



The first term gives the kinetic energy, and the second term gives the potential energy: this is consistent 
with our interpretation above, but now it is more general. Similarly, in Eq. ( 106| ) the "total momentum" 
consists of a mechanical term, and a potential term: 



p = (3^mc + qA. 



(108) 



The vector potential A contributes to the total momentum of the particle, in the same way that the 
scalar potential contributes to the total energy of the particle. Gauge invaiiance allows us to find new 
potentials that leave the fields (and hence the dynamics) of the system unchanged. Since the fields are 
obtained by taking derivatives of the potentials, this suggests that only changes in potentials between 
different positions and times are significant for the dynamics of charged particles. This in turn implies 
that only changes in (total) energy and (total) momentum are significant for the dynamics. 



7 Generation of electromagnetic waves 

As an example of the practical application of the potentials in a physical system, let us consider the gen- 
eration of electromagnetic waves from an oscillating, infinitesimal electric dipole. Although idealised, 
such a system provides a building block for constructing more realistic sources of radiation (such as the 
half-wave antenna), and is therefore of real interest. An infinitesimal electric dipole oscillating at a single 
frequency is known as an Hertzian dipole. 

Consider two point-like particles located on the z axis, close to and on opposite sides of the origin. 
Suppose that electric charge flows between the particles, so that the charge on each particle oscillates, 
with the charge on one particle being: 

qi = +goe-*"*, (109) 

and the charge on the other particle being: 

q2 = -qoe~"^'. (110) 
The situation is illustrated in Fig. [7] The current at any point between the charges is: 

I=^z = -iLoqoe~''''z. (Ill) 
at 

In the limit that the distance between the charges approaches zero, the charge density vanishes; however, 
there remains a non-zero electric current at the origin, oscillating at frequency lo and with amplitude Iq, 
where: 

Io = -iujqo- (112) 
Since the current is located only at a single point in space (the origin), it is straightforward to 



perform the integral in Eq. (98 1, to find the vector potential at any point away from the origin: 

A{r,t) = f^m- z, (113) 
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Fig. 7: Hertzian dipole: the charges oscillate around the origin along the z axis with infinitesimal amplitude. The 
vector potential at any point is parallel to the z axis, and oscillates at the same frequency as the dipole, with a phase 
difference and amplitude depending on the distance from the origin. 



where: 



k 



UJ 



(114) 



and £ is the length of the current: strictly speaking, we take the limit £ — 0, but we do so, we increase 
the current amplitude Iq, so that the produce IqI remains non-zero and finite. 

Notice that, with Eq. ( |1 13 1, we have quickly found a relatively simple expression for the vector 
potential around an Hertzian dipole. From the vector potential we can find the magnetic field; and from 
the magnetic field we can find the electric field. By working with the potentials rather than with the 
fields, we have greatly simplified the finding of the solution in what might otherwise have been quite a 
complex problem. 

The magnetic field is given, as usual, by S = V x ^. It is convenient to work in spherical polar 
coordinates, in which case the curl is written as: 



V X ^ 



1 



r-^ sin ( 



r rt) r sm U (p 

d_ d_ d_ 

dr ae d<t> 

Ar rAf) r sin 9A^ 



(115) 



Evaluating the curl for the vector potential given by Eq. ( 113 1 we find: 



Br = 0, 
Be = 0, 



(116) 
(117) 

(118) 



The electric field can be obtained from V x B = (which follows from Maxwells equation ^ in 



free space). The result is: 

Ep 



-^-(/o£)(l + f 
47reo c \ kr J 



(119) 
(120) 
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E4> = 0. (121) 

Notice that the expressions for the fields are considerably more complicated than the expression for the 
vector potential, and would be difficult to obtain by directly solving Maxwell's equations. 

The expressions for the fields all involve a phase factor e'^'^'^'^*), with additional factors giving 
the detailed dependence of the phase and amplitude on distance and angle from the dipole. The phase 
factor e*('^^~'^*) means that the fields propagate as waves in the radial direction, with frequency cu (equal 
to the frequency of the dipole), and wavelength A given by: 

A=^ = 5^. (122) 

k LO 

If we make some approximations, we can simplify the expressions for the fields. In fact, we can 
identify two different regimes. The near field regime is defined by the condition kr <^l.\n this case, the 
fields are observed at a distance from the dipole much less than the wavelength of the radiation emitted 
by the dipole. The dominant field components are then: 

« ^(V)sin^ (123) 

1 2i pi(kr-ujt) 

1 A „i{kr-ujt) 

Ee ^ (Jo^)-sing , „ . (125) 

47r£o c kr-^ 

The far field regime is defined by the condition A;r 2> 1. In this regime, the fields are observed 
at distances from the dipole that are large compared to the wavelength of the radiation emitted by the 
dipole. The dominant field components are: 

Un Jikr-ut) 

Ba, « -iV^(V) fcsin^ , (126) 

47r r 

1 k pi{kr-u>t) 

Ee -i- (Ioe)-sme . (127) 

Attsq c r 

The following features of the fields in this regime are worth noting: 

- The electric and magnetic field components are perpendicular to each other, and to the (radial) 
direction in which the wave is propagating. 

- At any position and time, the electric and magnetic fields are in phase with each other. 

- The ratio between the magnitudes of the fields at any given position is \E0\/\Bfp\ ^ c. 

These are all properties associated with plane waves in free space. Furthermore, the amplitudes of the 
fields falls off as 1/r: at sufficiently large distance from the oscillating dipole, the ampUtudes decrease 
slowly with increasing distance. At a large distance from an oscillating dipole, the electromagnetic waves 
produced by the dipole make a good approximation to plane waves in free space. 

It is also worth noting the dependence of the field amplitudes on the polar angle 9: the amplitudes 
vanish for ^ = and 6 = tt, i.e. in the direction of oscillation of the charges in the dipole. However, the 
ampUtudes reach a maximum for 9 = 7r/2, i.e. in a plane through the dipole, and perpendicular to the 
direction of oscillation of the dipole. 

We have seen that electromagnetic waves carry energy. This suggests that the Hertzian dipole 
radiates energy, and that some energy "input" will be required to maintain the amplitude of oscillation of 
the dipole. This is indeed the case. Let us calculate the rate at which the dipole will radiate energy. As 
usual, we use the expression: 

S = ExH, (128) 
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Fig. 8: Distribution of radiation power from an Hertzian dipole. The current in the dipole is oriented along the z 
axis. The distance of a point on the curve from the origin indicates the relative power density in the direction from 
the origin to the point on the curve. 



where the Poynting vector S gives the amount of energy in an electromagnetic field crossing unit area 
(perpendicular to S) per unit time. Before taking the vector product, we need to take the real parts of the 
expressions for the fields: 

= ^iIoi)k ^— ^ + sm{kr -cot)], (129) 

Att r \ kr J 

Er = -^hmU.o.{kr-u:t)- ''^^\-^'A , (130) 
47reo c \ kr J 

^ 1 „.k sinO ( sm( kr — u)t) cos( kr — ujt) ,\ 

= 4^'^°^'c — (- k^r^ + kr (»» 

The full expression for the Poynting vector will clearly be rather complicated; but if we take the average 
over time (or position), then we find that most terms vanish, and we are left with: 

As expected, the radiation is directional, with most of the power emitted in a plane through the dipole, 
and perpendicular to its direction of oscillation; no power is emitted in the direction in which the dipole 
oscillates. The power distribution is illustrated in Fig.[8] The power per unit area falls off with the square 
of the distance from the dipole. This is expected, from conservation of energy. 

The total power emitted by the dipole is found by integrating the power per unit area given by 



Eq. (132 1 over a surface enclosing the dipole. For simplicity, let us take a sphere of radius r. Then, the 



total (time averaged) power emitted by the dipole is 

I-2TT 



(P) = / \{S)\r'^ sine dOdcj). (133) 
Je=o J<j>=o 

Using the result: 



sin^ede=^, (134) 



we find: 



(V)!^! = M^. (135) 

Notice that, for a given amplitude of oscillation, the total power radiated varies as the square of the 
frequency of the oscillation. The consequences of this fact are famihar in an everyday observation. Gas 
molecules in the Earth's atmosphere behave as small oscillating dipoles when the electric charges within 
them respond to the electric field in the sunlight passing through the atmosphere. The dipoles re-radiate 
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surface current density Js 



Fig. 9: (a) Left: "Pill box" surface for derivation of the boundary conditions on the normal component of the 
magnetic flux density at the interface between two media, (b) Right: Geometry for derivation of the boundary 
conditions on the tangential component of the magnetic intensity at the interface between two media. 



the energy they absorb, a phenomenon known as Rayleigh scattering. The energy from the oscillating 
dipoles is radiated over a range of directions; after many scattering "events", it appears to an observer 
that the light comes from all directions, not just the direction of the original source (the sun). Eq. (135 1 
tells us that shorter wavelength (higher frequency) light is scattered much more strongly than longer 
wavelength (lower frequency) light. Thus, the sky appears blue. 



8 Boundary conditions 

Gauss' theorem and Stokes' theorem can be applied to Maxwell's equations to derive constraints on the 
behaviour of electromagnetic fields at boundaries between different materials. For RF systems in particle 
accelerators, the boundary conditions at the surfaces of highly-conductive materials are of particular 
significance. 



8.1 General boundary conditions 

Consider first a short cylinder or "pill box" that crosses the boundary between two media, with the flat 
ends of the cylinder parallel to the boundary, see Fig.|9](a). Applying Gauss' theorem to Maxwell's 
equation Q gives: 

/ V -BdV = I B-dS = 0, 

Jv JdV 

where the boundary dV encloses the volume V within the cyhnder. If we take the hmit where the length 
of the cylinder {2h - see Fig. |9](a)) approaches zero, then the only contributions to the surface integral 
come from the flat ends; if these have infinitesimal area dS, then since the orientations of these surfaces 
are in opposite directions on opposite sides of the boundary, and parallel to the normal component of the 
magnetic field, we find: 

- Bi^ dS + B2± dS = 0, 

where Bi± and B2± are the normal components of the magnetic flux density on either side of the bound- 
ary. Hence: 

Bi± = B2^. (136) 
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In other words, the normal component of the magnetic flux density is continuous across a boundary. 
Applying the same argument, but starting from Maxwell's equation ([T]l, we find: 

D21. - D,^ = ps, (137) 

where D± is the normal component of the electric displacement, and ps is the surface charge density 
(i.e. the charge per unit area, existing purely on the boundary). 

A third boundary condition, this time on the component of the magnetic field parallel to a bound- 
ary, can be obtained by applying Stokes' theorem to Maxwell's equation Q. In particular, we consider a 
surface S bounded by a loop dS that crosses the boundary of the material, see Fig.|9](b). If we integrate 
both sides of Eq. Q over that surface, and apply Stokes' theorem (|7]l, we find: 

V xH-dS= I H- dl= I J-dS+^ I D-dS, (138) 

S JdS Js Js 

where / is the total current flowing through the surface S. Now, let the surface S take the form of a thin 
strip, with the short ends perpendicular to the boundary, and the long ends parallel to the boundary. In the 
limit that the length of the short ends goes to zero, the area of S goes to zero: the electric displacement 
integrated over S becomes zero. In principle, there may be some "surfa ce cu rrent", with density (i.e. 



current per unit length) J^: this contribution to the right hand side of Eq. ( 138 1 remains non-zero in the 
limit that the lengths of the short sides of the loop go to zero. In particular, note that we are interested 
in the component of Jg that is perpendicular to the component of H parallel to the surface. We denote 
this component of the surface current density Js±. Then, we find from Eq. ( 138| ) (taking the limit of zero 
length for the short sides of the loop): 

H2ii-H^ = -J,^, (139) 

where is the component of the magnetic intensity parallel to the boundary at a point on one side of 
the boundary, and H2\\ is the component of the magnetic intensity parallel to the boundary at a nearby 
point on the other side of the boundar y. 

A final boundary condition can be obtained using the same argument that led to Eq. ( |139 l, but 
starting from Maxwell's equation Q. The result is: 

^2||=£^i|i, (140) 

that is, the tangential component of the electric field E is continuous across any boundary. 



8.2 Electromagnetic waves on boundaries 



The boundary conditions (136 1, (1371, (139 1, and (140 1 must be satisfied for the fields in an electro- 
magnetic wave incident on the boundary between two media. This requirement leads to the familiar 
phenomena of reflection and refraction: the laws of reflection and refraction, and the amplitudes of the 
reflected and refracted waves can be derived from the boundary conditions, as we shall now show. 

Consider a plane boundary between two media (see Fig.[TO)). We choose the coordinate system so 
that the boundary lies in the x-y plane, with the z axis pointing from medium 1 into medium 2. We write 
a general expression for the electric field in a plane wave incident on the boundary from medium 1 : 



Er 



(141) 



In order to satisfy the boundary conditions, there must be a wave present on the far side of the boundary, 
i.e. a transmitted wave in medium 2: 



Ex = Eoxe 



i{kj'-r—ujj't) 



(142) 
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Fig. 10: Incident, reflected, and transmitted waves on a boundary between two media. 



Let us assume that there is also an additional (reflected) wave in medium 1, i.e. on the incident side of 
the boundary. It will turn out that such a wave will be required by the boundary conditions. The electric 
field in this wave can be written: 

Er = EoRe'^''^-'^^'^''*\ (143) 

First of all, the boundary conditions must be satisfied at all times. This is only possible if all waves 
are oscillating with the same frequency: 

CJ/ = LOT = = w. (144) 

Also, the boundary conditions must be satisfied for all points on the boundary. This is only possible if 
the phases of all the waves vary in the same way across the boundary. Therefore, if p is any point on the 
boundary: 

ki ■ p = kr ■ p = kR- p. (145) 

Let us further specify our coordinate system so that kj lies in the x-z plane, i.e. the y component of kj 
is zero. Then, if we choose p to lie on the y axis, we see from Eq. ( 145 1 that: 

kry = kRy = kiy = 0. (146) 

Therefore, the transmitted and reflected waves also lie in the x-z plane. 

Now let us choose p to lie on the x axis. Then, again using Eq. ( 145 1, we find that: 

kTx = kRx = kj^. (147) 

If we define the angle Oj as the angle between kj and the z axis (the normal to the boundary), and 
similarly for 9t and 9r, then Eq. ( |147| ) can be expressed: 

kx sin 9t = kR sin Or = kj sin Oj. (148) 

Since the incident and reflected waves are travelling in the same medium, and have the same frequency, 
the magnitudes of their wave vectors must be the same, i.e. kR = kj. Therefore, we have the law of 
reflection: 

sin6lR = sin6l/. (149) 
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Fig. 11: Electric and magnetic fields in the incident, reflected, and transmitted waves on a boundary between two 
media. Left: The incident wave is N polarised, i.e. with the electric field normal to the plane of incidence. Right: 
The incident wave is P polarised, i.e. with the electric field parallel to the plane of incidence. 



The angles of the transmitted and incident waves must be related by: 

^ = 1^ = ^, (150) 
sm Ot ki V2 

where vi and V2 are the phase velocities in the media 1 and 2 respectively, and we have used the dispersion 
relation v = Lo/k. If we define the refractive index n of a medium as the ratio between the speed of light 
in vacuum to the speed of light in the medium: 

n=-, (151) 



then Eq. ( 150 1 can be expressed: 

^ = ^. (152) 

This is the familiar form of the law of refraction, Snell's law. 

So far, we have derived expressions for the relative directions of the incident, reflected, and trans- 
mitted waves. To do this, we have only used the fact that boundary conditions on the fields in the wave 
exits. Now, we shall derive expressions for the relative amplitudes of the waves: for this, we shall need 
to apply the boundary conditions themselves. 

It turns out that there are different relationships between the amplitudes of the waves, depending 
on the orientation of the electric field with respect to the plane of incidence (that is, the plane defined 
by the normal to the boundary and the wave vector of the incident wave). Let us first consider the case 
that the electric field is normal to the plane of incidence, i.e. 'W polarisation", see Fig.[TTJ left. Then, 
the electric field must be tangential to the boundary. Using the boundary condition ( 140[ ), the tangential 
component of the electric field is continuous across the boundary, and so: 

Eqi + Eon = EoT- (153) 



Using the boundary condition ( 139 1, the tangential component of the magnetic intensity H is also con- 
tinuous across the boundary. However, because the magnetic field in a plane wave is perpendicular to 
the electric field, the magnetic intensity in each wave must lie in the plane of incidence, and at an angle 
to the boundary. Taking the directions of the wave vectors into account: 

Hqj cos 6 1 — Hqji cos 9 1 = Hqt cos 6t- (154) 
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Using the definition of the impedance Z of a medium as the ratio between the ampHtude of the electric 
field and the amplitude of the magnetic intensity: 



Z 



El 

Ho' 



we can solve Eqs. (|153|) and (|154|) to give: 



/ Eor \ 
Eqi Jn 



Z2 cos 9 J — Z\ cos Bt 
Zi cos Bi + Z\ cos Qt ' 

2Z2 cos Qi 
Zi cos di + Z\ cos 



(155) 

(156) 
(157) 



Following a similar procedure for the case that the electric field is oriented so that it is parallel to 



the plane of incidence ("P polarisation", see Fig. 1 1 right), we find: 

Z2 cos Qt — Z\ cos Bi 
Z2 cos Ot + Z\ cos Qi ' 

2Z2 cos Qi 
Z2 cos Qt + Z\ cos 0/ 



/ Eqr \ 
Eqi Jp 



(158) 
(159) 



Equations ( 156 1, ( 157 1, ( 158 1 and ( 159 1 are known as Fresnel's equations: they give the amplitudes 
of the reflected and transmitted waves relative to the amplitude of the incident wave, in terms of the 
properties of the media (specifically, the impedance) on either side of the boundary, and the angle of 
the incident wave. Many important phenomena, including total internal reflection, and polarisation by 
reflection, follow from Fresnel's equations. However, we shall focus on the consequences for a wave 
incident on a good conductor. 

First, note that for a dielectric with permittivity e and permeability fi, the impedance is given by: 



Z 



(160) 



This follows from Eq. ( 155 1, using the constitutive relation B = fiH, and the relation between the field 
amplitudes in an electromagnetic wave Eq/Bq = v, where the phase velocity v = l/^/JIs. 

Now let us consider what happens when a wave is incident on the surface of a conductor. Using 
Maxwell's equation (|4]), the impedance can be written (in general) for a plane wave: 



Z 



Eo 
Ho 



(161) 



In a good conductor (for which the conductivity a » Loe), the wave vector is complex; and this means 
that the impedance will also be complex. This implies a phase difference between the electric and 
magnetic fields, which does indeed occur in a conductor. In the context of Fresnel's equations, complex 
impedances will describe the phase relationships between the incident, reflected, and transmitted waves. 



From Eq. (52 1, the wave vector in a good conductor is given (approximately) by: 



Therefore, we can write: 



Z 




(1 




(162) 



(163) 
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Consider a wave incident on a good conductor from a dielectric. If the permittivity and permeability 
of the conductor are similar to those in the dielectric, then, since a ^ uje (by definition, for a good 
conductor), the impedance of the conductor will be much less than the impedance of the dielectric. 
Fresnel's equations become (for both and P polarisation): 

« -1, - 0. (164) 

There is (almost) perfect reflection of the wave (with a change of phase); and very little of the wave 
penetrates into the conductor. 

At optical frequencies and below, most metals are good conductors. In practice, as we expect from 
the above discussion, most metals have highly reflective surfaces. This is of considerable importance for 
RF systems in particle accelerators, as we shall see when we consider cavities and waveguides, in the 
following sections. 



8.3 Fields on the boundary of an ideal conductor 

We have seen that a good conductor will reflect most of the energy in a wave incident on its surface. 
We shall define an ideal conductor as a material that reflects all the energy in an electromagnetic wave 
incident on its surfaced In that case, the fields at any point inside the ideal conductor will be zero at all 



times. From the boundary conditions ( 136 1 and ( 140 1, this implies that, at the surface of the conductor: 



Si 



0, 



^11 = 0. 



(165) 



That is, the normal component of the magnetic field, and the tangential component of the electric field 
must vanish at the boundary. These conditions impose strict constraints on the patterns of electromag- 
netic field that can persist in RF cavities, or that can propagate along waveguides. 



The remaining boundary conditions, ( 137 1 and ( 139 1, allow for discontinuities in the normal com- 
ponent of the electric field, and the tangential component of the magnetic field, depending on the presence 
of surface charge and surface current. In an ideal conductor, both surface charge and surface current can 
be present: this allows the field to take non-zero values at the boundary of (and within a cavity enclosed 
by) an ideal conductor. 



9 Fields in cavities 

In the previous section, we saw that most of the energy in an electromagnetic wave is reflected from 
the surface of a good conductor. This provides the possibility of storing electromagnetic energy in the 
form of standing waves in a cavity; the situation will be analogous to a standing mechanical wave on, 
say, a violin string. We also saw in the previous section that there are constraints on the fields on the 
surface of a good conductor: in particular, at the surface of an ideal conductor, the normal component 
of the magnetic field and the tangential component of the electric field must both vanish. As a result, 
the possible field patterns (and frequencies) of the standing waves that can persist within the cavity are 
determined by the shape of the cavity. This is one of the most important practical aspects for RF cavities 
in particle accelerators. Usually, the energy stored in a cavity is needed to manipulate a charged particle 
beam in a particular way (for example, to accelerate or deflect the beam). The effect on the beam is 
determined by the field pattern. Therefore, it is important to design the shape of the cavity, so that the 
fields in the cavity interact with the beam in the desired way; and that undesirable interactions (which 
always occur to some extent) are minimised. The relationship between the shape of the cavity and the 
different field patterns (or modes) that can persist within the cavity will be the main topic of the present 
section. Other practical issues (for example, how the electromagnetic waves enter the cavity) are beyond 
our scope. 

'it is tempting to identify superconductors with ideal conductors; however, superconductors are rather complicated materi- 
als, that show sometimes surprising behaviour not always consistent with our definition of an ideal conductor. 
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Fig. 12: Rectangular cavity. 



9.1 Modes in a rectangular cavity 

We consider first a rectangular cavity with perfectly conducting walls, containing a perfect vacuum (see 



Fig. 12 1. The wave equation for the electric field inside the cavity is: 



where c is the speed of light in a vacuum. There is a similar equation for the magnetic field B. We 
look for solutions to the wave equations for E and B that also satisfy Maxwell's equations, and also 
satisfy the boundary conditions for the fields at the walls of the cavity. If the walls of the cavity are ideal 
conductors, then the boundary conditions are: 

Ei\ = 0, (167) 
B_i_ = 0, (168) 

where E'n is the component of the electric field tangential to the wall, and B± is the component of the 
magnetic field normal to the wall. 

Free-space plane wave solutions will not satisfy the boundary conditions. However, we can look 
for solutions of the form: 

E{x,y,z,t) = Er^e-'^^^ (169) 

where Ef^ = Ef{x, y, z) is a vector function of position (independent of time). Substituting into the wave 
equation, we find that the spatial dependence satisfies: 

V'^Ef^+'^Ef^ = 0. (170) 

The full solution can be derived using the method of separation of variables (in fact, we have begun the 
process by separating the time from the spatial variables). However, it is sufficient to quote the result, 
and verify the solution simply by substitution into the wave equation. The components of the electric 
field in the rectangular cavity are given by: 

Ex = Exo cos kxX sin kyU sin kzZ e~^^^, (171) 
Ey = Eyo sin kxX cos kyU sin kzze~^'^^, (172) 



Ez = Ezo sin kxX sin kyU cos kzZ e (173) 



29 



Fig. 13: Boundary conditions in a rectangular cavity. The field component Ey is parallel to the walls of the cavity 
(and must therefore vanish) at a; = and x — a^- 



To satisfy the wave equation, we require: 

kj, -\- hy -\- = 

Maxwell's equation V ■ E = 0, imposes a constraint on the components of the wave vector and the 
amplitudes of the field components: 

kxExo + kyEyo + k;,E^o = 0. (175) 

We also need to satisfy the boundary conditions, in particular that the tangential component of the 
electric field vanishes at the walls of the cavity. This imposes additional constraints on k^, ky and k^. 
Consider: 

Ey = Eyo sin kxX cos kyt/ sin kzZ e~*'^*. (176) 



The boundary conditions require that Ey = at x = and x = ax, for all y, z, and t (see Fig. 13 1. These 
conditions are satisfied if kxUx = rrixTT, where nix is an integer. To satisfy all the boundary conditions, 
we require: 

Kx , Ivy , Kz , ( i ' ' ) 

QjX CLy dz 

where nix, "my and are integers. Note that these integers play a large part in determining the shape of 
the electric field (though even when these numbers are specified, there is still some freedom in choosing 
the relative amplitudes of the field components). The quantities nix, iT^y, and niz, are called the mode 
numbers. The frequency of oscillation is determined completely by the mode numbers, for a given size 
and shape of cavity. 

The magnetic field can be obtained from the electric field, using Maxwell's equation: 

dB 

VxE = - — . (178) 

This gives: 

Ezoky) sin kxX cos ky-y cos kzZ e~*'^*, (179) 
Exokz) coskxX sinkyy cos kzZ e~'^'^^, (180) 

EyQkx) coskxX coskyy svukzZ e'"^^^ . (181) 

It is left as an exercise for the reader to show that these fields satisfy the boundary condition on the 
magnetic field at the walls of the cavity, and also satisfy the remaining Maxwell's equations: 



Bx = 


1 

— {Eyokz 


By = 


-{Ezokx 


OJ 


Bz = 


i 

— [Exoky 



1 f)p 

V-B = 0, and V x B = ^ — . (182) 

ot 
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0.5 1 1.5 2 2.5 3 

to/jtC 

Fig. 14: Mode spectra in rectangular cavities. Top: all side lengths equal. Middle: two side lengths equal. Bottom: 
all side lengths different. Note that we show all modes, including those with two (or three) mode numbers equal 
to zero, even though such modes will have zero amplitude. 



Note that the frequency of oscillation of the wave in the cavity is determined by the mode numbers 

nix, iTT'D and nriz'. 

u = J(^)\(^)\(^)\ (183) 

For a cubic cavity (a^ = ay = a^), there will be a degree of degeneracy, i.e. there will generally be 
several different sets of mode numbers leading to different field patterns, but all with the same frequency 



of oscillation. The degeneracy can be broken by making the side lengths different: see Fig. 14 



Some examples of field patterns in different modes of a rectangular cavity are shown in Fig. 15 



9.2 Quality factor 

Note that the standing wave solution represents an oscillation that will continue indefinitely: there is 
no mechanism for dissipating the energy. This is because we have assumed that the walls of the cavity 
are made from an ideal conductor, and the energy incident upon a wall is completely reflected. These 
assumptions are implicit in the boundary conditions we have imposed, that the tangential component of 
the electric field and the normal component of the magnetic field vanish at the boundary. 

In practice, the walls of the cavity will not be perfectly conducting, and the boundary conditions 
will vary slightly from those we have assumed. The electric and (oscillating) magnetic fields on the walls 
will induce currents, which will dissipate the energy. The rate of energy dissipation is usually quantified 
by the "quality factor", Q. The equation of motion for a damped harmonic oscillator: 



d u uj du 
~d£^ ^ Qdt 

has the solution: 



+ w\ = (184) 



u = uqc 2Q cos(a;'t — i^), (185) 

where: 

The quality factor Q is (for Q ^ I) the number of oscillations over which the energy in the oscillator 
(proportional to the square of the amplitude of u) falls by a factor 1/e. More precisely, the rate of energy 
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Fig. 15: Examples of modes in rectangular cavity. From top to bottom: (110), (111), (210), and (211). 
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dissipation (the dissipated power, P^) is given by: 

If a field is generated in a rectangular cavity corresponding to one of the modes we have calculated, 
the fields generating currents in the walls will be small, and the dissipation will be slow: such modes 
(with integer values of rux, niy and m^) will have a high quality factor, compared to other field patterns 
inside the cavity. A mode with a high quaUty factor is called a resonant mode. In a rectangular cavity, 
the modes corresponding to integer values of the mode numbers are resonant modes. 

9.3 Energy stored in a rectangular cavity 

It is often useful to know the energy stored in a cavity. For a rectangular cavity, it is relatively straight- 
forward to calculate the energy stored in a particular mode, given the amplitude of the fields. The energy 
density in an electric field is: 

Ue = ^D- E. (188) 
Therefore, the total energy stored in the electric field in a cavity is: 

£E = \eojE^dV, (189) 

where the volume integral extends over the entire volume of the cavity. In a resonant mode, we have: 

/ cos kxxdx= / sin kxxdx=-, (190) 
Jo Jo 2 

where kx = mxTr/ax, and nix is a non-zero integer. We have similar results for the y and z directions, 
so we find (for rUx, ruy and all non-zero integers): 

£e = ^eo{Elo + E^o + E^zo) cos' ^t. (191) 

The energy varies as the square of the field amplitude, and oscillates sinusoidally in time. 
Now let us calculate the energy in the magnetic field. The energy density is: 

Ub = Ib- H. (192) 

kl + kl + k'^^ = and Exokx + Eyoky + E^ok^, = 0, (193) 

we find, after some algebraic manipulation (and noting that the magnetic field is 90° out of phase with 
the electric field): 

£b = ^-^ {Elo + + eIq) sit? ojt. (194) 

As in the case of the electric field, the numerical factor is correct if the mode numbers rUx, my and ruz 
are non-zero integers. 

Finally, using 1/c? = fioSo, we have (for nix, my and non-zero integers): 

£e + £b = ^eo {Elo + E^yo + eIq). (195) 

The total energy in the cavity is constant over time, although the energy "oscillates" between the electric 
field and the magnetic field. 



Using: 
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The power flux in the electromagnetic field is given by the Poynting vector: 

S = ExH. (196) 

Since the electric and magnetic fields in the cavity are 90° out of phase (if the electric field varies as 
cosujt, then the magnetic field varies as sinc^t), averaging the Poynting vector over time at any point in 
the cavity gives zero: this is again consistent with conservation of energy. 



9.4 Shunt impedance 

We have seen that in practice, some of the energy stored in a cavity will be dissipated in the walls, and 



that the rate of energy dissipation for a given mode is measured by the quality factor, Q (Eq. ( 1 87 1) 



dt~ Q^- 

For a mode with a longitudinal electric field component Ezo = Vq/L (where L is the length of the 
cavity), we define the shunt impedance, Rg-. 



T/2 
2Pd 



(197) 



(Note that different definitions of the shunt impedance are used, depending on the context: some defini- 
tions omit the factor 1/2). 



Combining Eqs. ( 187 1 and ( 197 1, we see that 



Rs _ Vq Pel 



1^2 



(198) 



Q 2Pd io£ 2uj£ 

Consider a mode with = 0. Such modes have only transverse components of the magnetic 



field, and are called TM modes. Using equation ( 181 1, we see that the electric field in TM modes obeys: 

kyE^Q = k^EyO. (199) 



We also have, from ( 175 1: 



k^E^ + kyEyo + k,E,o = 0. (200) 

These relations allow us to write the energy stored in the cavity purely in terms of the mode numbers and 
the peak longitudinal electric field: 



Combining equations ( |198| ) and ( |201| ), we see that: 

U2 I 1.2 



Q 



^X + ky 



L2 



(201) 



(202) 



For a TM mode in a rectangular cavity, the quantity Rg/Q depends only on the length of the cavity 
and the mode numbers. It is independent of such factors as the material of the walls. In fact, this result 
generalises: for TM modes, Rg/Q depends only on the geometry of the cavity, and the mode numbers. 
This is of practical significance since, to optimise the design of a cavity for accelerating a beam, the goal 
is to maximise Rg/Q for the accelerating mode, and minimise this quantity for all other modes. Since 
Rg/Q is independent of such quantities as the material from which the cavity is made, it is possible while 
designing the cavity to focus on optimising the geometry to achieve the highest Rg/Q for the desired 
mode (and minimise Rg/Q) for the other modes. Properties of the material from which the cavity will 
be made can safely be neglected at this stage of the design process. 
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9.5 Cylindrical cavities 

Most cavities in accelerators are closer to a cylindrical than a rectangular geometry. It is worth consid- 
ering the solutions to Maxwell's equations, subject to the usual boundary conditions, for a cylinder with 
perfectly conducting walls. We can find the modes in just the same way as we did for a rectangular cav- 
ity: that is, we find solutions to the wave equations for the electric and magnetic fields using separation 
of variables; then we find the "allowed" solutions by imposing the boundary conditions. The algebra 
is more complicated this time, because we have to work in cylindrical polar coordinates. We will not 
go through the derivation in detail: the solutions for the fields can be checked by taking the appropriate 
derivatives. 

One set of modes (not the most general solution) we can write down is as follows: 

Er = -iBo^JJkrr) sinnd sink^ze-''^^ (203) 

Eg = -iBo^4X^rr) cosnO smk^.ze'''^^ (204) 

E, = 0, (205) 

k 

Br = Bo-r^Jnikrr) COS n0 COS kzze~"^\ (206) 

Bg = -Bo^JJkrr) sinnO COS k^ze'''^^ (207) 
k^r 

B^ = Bo J„ikrr) cos nO sin k^^ze'^''^ (208) 
Note that Jn{x) is a Bessel function of order n, and Jn{x) is the derivative of Jn{x). The Bessel functions 



(Fig. 16 1 are solutions of the differential equation: 



y" + - + fi-^')y = o. (209) 

X \ X'^ ) 

This equation appears when we separate variables in finding a solution to the wave equation in cylindrical 
polar coordinates. 

Because of the dependence of the fields on the azimuthal angle 6*, we require that n is an integer: 
n provides an azimuthal index in specifying a mode. 

From the boundary conditions. Eg and B^ must both vanish on the curved wall of the cavity, i.e. 
when r = a, where a is the radius of the cylinder. Therefore, we have a constraint on kr'. 

J^{kra) = 0, (210) 

or: ^ 

kr = —, (211) 

a 

where is the mth zero of the derivative of the nth order Bessel function. This equation is analogous 
to the conditions we had for the rectangular cavity, e.g. k^ = nix'K/ax- We can use the integer m as a 
radial index in specifying a mode. 

We also need to have B^ = Er = Eg = on the flat ends of the cavity. Assuming the flat ends of 
the cavity are at z = and z = L, these boundary conditions are satisfied if: 

in 

sinkzL = 0, therefore kz = —, (212) 
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Fig. 16: Bessel functions. The red, green, blue and purple lines show (respectively) the first, second, third and 
fourth order Bessel functions J„(a;). 



where ^ is an integer. £ provides a longitudinal index in specifying a mode. 
Also, we find that: 



V'^E = -{kl + kl)E, (213) 



so from the wave equation: 



we must have: 



~\~ k~ — . (215) 



Similar equations hold for the magnetic field, B. 

In the modes that we are considering, the longitudinal component of the electric field = 0, i.e. 
the electric field is purely transverse. These modes are known as "TE" modes. A specific mode of this 
type, with indices n (azimuthal), m (radial), and £ (longitudinal) is commonly written as TEnmi- The 
frequency of mode TEj-imi depends on the dimensions of the cavity, and is given by: 



c/WT^ = cJ + (^y. (216) 



<^nm£ 



The fields in the TEno mode in a cylindrical cavity are shown in Fig. 17 



TE modes are useful for giving a transverse deflection to a beam in an accelerator, but are not 
much use for providing acceleration. Fortunately, cylindrical cavities allow another set of modes that 
have non-zero longitudinal electric field: 

k 

Er = -Eoj^ Jnikrr) COS nO sin k^ze-"^\ (217) 

Eg = Eo^Jnikrr) sinnO smk,ze~''^\ (218) 

E^ = EoJn{krr) cos ne cos k;,ze~''^\ (219) 
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Fig. 17: TEiio mode in a cylindrical cavity. Left: Electric field. Right; Magnetic field. Top: 3 -dimensional view. 
Bottom: Cross-sectional view. 



Br = iEQ^—r-Jn{krr) sinnO cosk-^ze~''^\ (220) 
Be = iEo-^J'^ikrr) cosnO cosk^ze'''^^ (221) 

B^ = 0. (222) 

In these modes, the magnetic field is purely transverse (zero longitudinal component); therefore, they are 
referred to as "TM" modes. As before, for physical fields, n must be an integer. The boundary conditions 
on the fields give: 

kr = and k^ = (223) 

a L 

where pnm is the mth zero of the nth order Bessel function Jn{x). The frequency of a mode TM^m^ is 
given by: 

Unml = Cy^WTM = (P^f + . (224) 

The lowest frequency accelerating mode in a cylindrical cavity is the TMqio mode (n = 0, m = 1, 



£ = 0). The fields in the TMqio (see Fig. 18 1 mode are given by: 



Er = 0, (225) 
Ee = 0, (226) 



E, = EoJoiPoi-) e~'"N (227) 



Br = 0, (228) 
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Fig. 18: Fields in the TMgio mode in a cylindrical cavity. Left: Electric field. Right: Magnetic field. 
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iujt 



0. 



(229) 
(230) 



The frequency of the TMqio mode is determined by the radius of the cavity, and is independent of 
its length: 

woio=Poi-- (231) 
a 

(Note that poi ~ 2.40483.) However, to get the maximum acceleration from the cavity, the time taken 
for a particle to pass through the cavity should be one half of the RF period, i.e. vr/w. Therefore, for 
best efficiency, the length of the cavity should be L = vru/o;, where v is the velocity of the particle. For 
ultra-relativistic particles, the length of the cavity should be L = A/2, where A is the wavelength of an 
electromagnetic wave with angular frequency a; in free space. 



10 Waveguides 

Cavities are useful for storing energy in electromagnetic fields, but it is also necessary to transfer electro- 
magnetic energy between different locations, e.g. from an RF power source such as a klystron, to an RF 
cavity. Waveguides are generally used for carrying large amounts of energy (high power RF). For low 
power RF signals (e.g. for timing or control systems), transmission lines are generally used. Although 
the basic physics in waveguides and transmission lines is the same - both involve electromagnetic waves 
propagating through bounded regions - different formalisms are used for their analysis, depending on 
the geometry of the boundaries. In this section, we consider waveguides; we consider transmission lines 
in the following section. 

As was the case for cavities, the patterns of the fields in the resonant modes are determined by 
the geometry of the boundary. The resonant modes in a waveguide with circular cross-section will 
look somewhat different from the modes in a waveguide with rectangular cross-section. The physical 
principles and general methods are the same in both cases; however, the algebra is somewhat easier for 
a rectangular waveguide. Therefore, we shall consider here only waveguides with rectangular cross- 
section. 



10.1 Modes in a rectangular waveguide 

Consider an ideal conducting tube with rectangular cross-section, of width ax and height ay (Fig. 19 1. 
This is essentially a cavity resonator with length — )• oo. The electric field must solve the wave 
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Fig. 19: Rectangular waveguide. 



equation: 



together (as usual) with Maxwell's equations. By comparison with the rectangular cavity case, we expect 
to find standing waves in x and y, with plane wave solution in z. Therefore, we try a solution of the form: 

Ex = ExQCoskxX sinkyye'^^''''~'^*\ (233) 

Ey = Eyo sin kxx cos kyy e^^*^-^"'^*) , (234) 

E, = iE^osmkxX s\nkyye'^^'''-'^^\ (235) 
To satisfy the wave equation, we must have: 



1 d'^E 



kl + kl + kl = 'jj. (236) 

We must of course also satisfy Maxwell's equation: 

V • = 0, (237) 

at all X, y and z. This leads to the condition: 

kxExo + kyEyo + kzE;,o = 0, (238) 

which is the same as we found for a rectangular cavity. 

Now we apply the boundary conditions. In particular, the tangential component of the electric 
field must vanish at all boundaries. For example, we must have Ez = for all t and all z, at x = and 



2 



X = ax, and at y = and y = ay (see Fig. 19 1. Therefore, sin kxax = 0, and sin A;„a„ = 0, and hence: 



kx = — , m^. = 0,1,2... (239) 

ax 

m,,7r 

ky = my = 0,1, 2... (240) 
ay 

These conditions also ensure that Ex = at y = and at y = ay, and that Ey = at x = and at 

X = Or- 
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We can find the magnetic field in tlie waveguide from: 



V xE 



dB 

Ik' 



The result is: 



= -{kyE-,Q - kzEyo) sin k^x cos kyy e'^^'^ '^^^ , 
B„ = -{k^EM - k^E^o) cosk^x sin A;^?/ e*(^-^~'^*\ 



CO 



B, 



-{k^EyQ - kyE^) COS cos kyy e'^''"'' 



lUJ 



(241) 

(242) 
(243) 
(244) 



The conditions (239 1 and (240i: 



nirTT 



_ myn 



0,1,2. 



0,1,2. 



ensure that the boundary conditions on the magnetic field are satisfied. Also, the conditions (236) 



and (238 1: 



2,2 I 2,2 , -l2 _ ^ 



kxExo + kyEyo + kzEzo 



0, 



ensure that Maxwell's equation: 



V X B 



1 dE 



(245) 



is satisfied. 

In a rectangular cavity, the mode could be specified by a set of three numbers, nix, f^y, and m^, 
that determined the wave numbers and the frequency of oscillation of the field. Because nix, my and niz 
had to be integers, there was a discrete spectrum of allowed frequencies of oscillation in a cavity. In a 
waveguide, we only need two integer mode numbers, and niy, corresponding to the wave numbers 
in the transverse dimensions. The third wave number, in the longitudinal direction, is allowed to take 
values over a continuous range: this is because of the absence of boundaries in the longitudinal direction. 
Therefore, there is a continuous range of frequencies also allowed. 

When we discussed cylindrical cavities, we wrote down two sets of modes for the fields: TE modes 
{Ez = 0), and TM modes {Bz = 0). We can also define TE and TM modes in rectangular cavities, and 
rectangular waveguides. For a TE mode in a waveguide, the longitudinal component of the electric field 
is zero at all positions and times, which implies that: 



EzG = 0. 



(246) 



From (238 1, it follows that in a TE mode, the wave numbers and amplitudes must be related by: 

kxExo + kyEyQ = 0. (247) 



Similarly, for a TM mode, we have (using 244 1: 

Bzo = 0, kxEyo — kyExo = 0. 



(248) 
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10.2 Cut-off frequency in a waveguide 

We have seen that in a waveguide, specifying the mode does not uniquely specify the frequency. This is 
a consequence of the fact that there are no boundaries in the longitudinal direction: there are only two 
mode numbers required to specify the transverse dependence of the field; the longitudinal dependence is 
not constrained by boundaries. The frequency of oscillation of fields in the waveguide depends on the 
longitudinal variation of the field, as well as on the transverse variation. However, if we want a wave 
to propagate along the waveguide, there is a minimum frequency of oscillation that is needed to achieve 
this, as we shall now show. 



Let us write equation (236 1 in the form: 



^2 

kz — k^ — ky. (249) 



If the wave is to propagate with no attentuation (i.e. propagate without any loss in amplitude), then kz 
must be real. We therefore require that: 



OJ 



2 



>kl + kl. (250) 



c 



2 " '"x > '"y 



This can be written: 



We define the cut-off frequency locq: 




(251) 




Wco = ^cW — + . (252) 



Above the cut-off frequency, the wave will propagate without loss of amplitude (in a waveguide with 
ideal conducting walls). At the cut-off frequency, the group velocity of the wave is zero: the wave 
is effectively a standing wave. Below the cut-off frequency, oscillating fields can still exist within the 
waveguide, but they take the form of attenuated (evanescent) waves, rather than propagating waves: the 
longitudinal wave number is purely imaginary. 

Note that at least one of the mode numbers mx and my must be non-zero, for a field to be present 
at all. Therefore, the minimum cut-off frequency of all the possible modes is given by: 

TT 

wco = -c, (253) 

a 

where c is the speed of light, and a is the side length of the longest side of the waveguide cross-section. 



10.3 Phase velocity, dispersion, and group velocity 

Consider an electromagnetic wave that propagates down a waveguide, with longitudinal wave number kz 
and frequency uj. The phase velocity Vp is given by: 




(254) 



This is the speed at which one would have to travel along the waveguide to stay at a constant phase of 
the electromagnetic wave (i.e. kzZ — Lot = constant). However, since kx, ky and kz are all real, it follows 
that the phase velocity is greater than the speed of light: 

Vp > c. (255) 
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Fig. 20: Dispersion curve (in red) for waves in a waveguide. The broken black line shows phase velocity uj/kz ~ c. 



The dispersion curve (Fig. 20 1 shows how the frequency uj varies with the longitudinal wave number kz- 

Although the phase velocity of a wave in a waveguide is greater than the speed of light, this does 
not violate special relativity, since energy and information signals travel down the waveguide at the group 
velocity, rather than the phase velocity. The group velocity Vg is given by: 



V, 



so we have: 



kz kz 2 

—c = — c , 

U2 \ U2 ih.2 00 



Vg < C. 



Note that for a rectangular waveguide, the phase and group velocities are related by: 

VpVg = (?. 



(256) 



(257) 



(258) 



Using equation (236 1: 



2,2 I 2,2 , z,2 _ ^ 



and the boundary conditions ( 239 1 and ( 240 1 



rriTiT 



and ky = — —, 



we can write the group velocity (256 1 for a given mode {mx, my) in terms of the frequency w: 



7r2c2 / ml "T.2 
w2 \al al 



(259) 



Energy in the wave propagates along the waveguide only for uj > Wco (where Wco is the cut-off fre- 
quency). Note the limiting behaviour of the group velocity, for oj > ujco- 



lim Vn 



(260) 
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lim Vq = 0. (261) 



The variation of group velocity with frequency is illustrated in Fig. 21 



10.4 Energy density and energy flow in a waveguide 

The main purpose of a waveguide is to carry electromagnetic energy from one place to another: therefore, 
it is of interest to know the energy density and the energy flow in a waveguide, for a given mode and field 
amplitude. 

The energy density in the waveguide can be calculated from: 

U=\eoE^ + \ — . (262) 



Using the expressions for the fields (233 1 - (235 1, and assuming that both kx and ky are non-zero, we 



see that the time-average energy per unit length in the electric field is: 

{Ue) dx dy = ^eoA {eI^ + E^^ + ^'o) > (263) 

where A is the cross-sectional area of the waveguide. 

The fields in a waveguide are just a superposition of plane waves in free space. We know that in 
such waves, the energy is shared equally between the electric and magnetic fields. Using the expressions 



for the magnetic field (242 1 - (244 1, we can indeed show that the magnetic energy is equal (on average) 
to the electric energy - though the algebra is rather complicated. The total time-average energy per unit 
length in the waveguide is then given by: 

(U) dx dy = '^eoA [eI, + EI^ + El,) . (264) 

The energy flow along the waveguide is given by the Poynting vector, S = E x H. Taking the 
real parts of the field components, we find for the time-average transverse components: 

{Sx) = {Sy) = 0. (265) 

The longitudinal component is non-zero; taking the time average and integrating over the cross section 
gives the time average energy flow along the waveguide 



{S,)dxdy = ^A{E^,o + El, + El,) ^ . (266) 
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Fig. 22: Inductance in a transmission line. 



Using equation (256 1 for the group velocity Vg, the time-average energy flow along the waveguide can 



be written in the form: 



(5,) dx dy = \eoA{El^ + E^^ + El^)vg. (267) 



8 

Thus, the energy density per unit length and the rate of energy flow are related by: 



{S^) dxdy = vg j j {U) dx dy. (268) 

This is consistent with our interpretation of the group velocity as the speed at which energy is carried by 
a wave. 

The boundary conditions that we have applied implicitly assume ideal conducting walls. In that 
case, there will be no energy losses to the walls, and electromagnetic waves can propagate indefinitely 
along the waveguide without any reduction in amplitude. In practice, there will be some attenuation, 
because of the finite conductivity of the walls of the wavguide. However, using a metal with high 
conductivity (such as copper or aluminium), attenuation lengths can be quite long, of order many metres. 
This makes waveguides very suitable for high-power RF applications. 



11 Transmission lines 

Transmission lines are used (as are waveguides) to guide electromagnetic waves from one place to an- 
other. A coaxial cable (used, for example, to connect a radio or television to an aerial) is an example of 
a transmission line. Transmission lines may be less bulky and less expensive than waveguides; but they 
generally have higher losses, so are more appropriate for carrying low-power signals over short distances. 



11.1 LC model of a transmission line 



Consider an infinitely long, parallel wire with zero resistance. In general, the wire will have some 
inductance per unit length, L, which means that when an alternating current / flows in the wire, there 



will be a potential difference between different points along the wire (Fig. 22 ). If V is the potential at 
some point along the wire with respect to earth, then the potential difference between two points along 
the wire is given by: 

'>^' (269) 



Ay 



dx 



6x 



-L6x 



dt 



In general, as well as the inductance, there will also be some capacitance per unit length, C, 



between the wire and earth (Fig. 23 ). This means that the current in the wire can vary with position: 



dx 



Sx 



-C6x 



dV 
~dt' 



(270) 
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Fig. 23: Capacitance in a transmission line. 



The "complete" model, including both inductance and capacitance, is referred to as the LC model for a 
transmission line. 

We wish to solve the equations for the current and voltage in the transmission line, as functions of 

(271) 



position and time. Let us take Equations (269 1 and (270 1 above: 

dl 



dV 
dx 



-L 



m 

dx 



-C 



dV 



(272) 



These equations are analogous to Maxwell's equations: they show how changes (with respect to time and 
position) in current and voltage are related to each other. They take the form of two, coupled first-order 
differential equations. The fact that they are coupled means that the solution is not obvious; however, 
we can take the same approach that we did with Maxwell's equations in free space. By taking additional 
derivatives, we can arrive at two uncoupled second-order differential equations. The solution to the 



uncoupled equations is more apparent than the solution to the coupled equations. Differentiate (271 
with respect to t: 



and ( 272 1 with respect to x\ 



Hence: 



dxdt 

c^l _ 
dx"^ 

dx"^ 



-L 



-c 



LC 



dtdx 



Similarly (by differentiating (271 1 with respect to x and (272 1 with respect to t), we find: 



dx"^ 



LC 



• 



(273) 



(274) 



(275) 



(276) 



Equations (275 1 and (276 1 are wave equations for the current in the wire, and the voltage between 
the wire and earth. The waves travel with speed v, given by: 



1 



The solutions to the wave equations may be written: 



(277) 



(278) 
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Fig. 24: Termination of a transmission line with an impedance R. 



Joe 



i{kx—uit) 



(279) 



where the phase velocity is: 



k 



(280) 



Note that the inductance per unit length L and the capacitance per unit length C are real and positive. 
Therefore, if the frequency a; is real, the wave number k will also be real: this implies that waves prop- 
agate along the transmission line with constant amplitude. This result is expected, given our assumption 
about the line having zero resistance. 



The solutions must also satisfy the first-order equations (269 1 and (270 1. Substituting the above 
solutions into these equations, we find: 



kVo = ojLIq, 
klo = ujCVq. 



(281) 
(282) 



Hence: 

The ratio of the voltage to the current is called the characteristic impedance, Z, of the transmission line. 
Z is measured in ohms, Q. Note that, since L and C are real and positive, the impedance is a real number: 
this means that the voltage and current are in phase. The characteristic impedance of a transmission line 
is analogous to the impedance of a medium for electromagnetic waves: the impedance of a transmission 
line gives the ratio of the voltage amplitude to the current amplitude; the impedance of a medium for 
electromagnetic waves gives the ratio of the electric field amplitude to the magnetic field amplitude. 



11.2 Impedance matching 

So far, we have assumed that the transmission line has infinite length. Obviously, this cannot be achieved 
in practice. We can terminate the transmission line using a "load" with impedance Zl that dissipates 
the energy in the wave while maintaining the same ratio of voltage to current as exists all along the 



transmission line - see Fig. 24 In that case, our above analysis for the infinite line will remain valid for 
the finite line, and we say that the impedances of the line and the load are properly matched. 

What happens if the impedance of the load, Zl, is not properly matched to the characteristic 
impedance of the transmission line, Zl In that case, we need to consider a solution consisting of a 
superposition of waves travelling in opposite directions: 

V = Y^e'^kx-uji) ^ j^y^^i(-kx-ujt) ^ (-284) 
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Fig. 25: A "lossy" transmission line. 



The corresponding current is given by: 



J = ^0 ^i(kx-ujt) _ j^yO_^i{-kx-u]t) ^ (285) 

Z Z 



Note the minus sign in the second term in the expression for the current: this comes from equations (271 



and (272 1. Let us take the end of the transmission line, where the load is located, to be at x = 0. At this 



position, we have: 

V = Vbe"'""* (1 + i^) , (286) 



/ = Y^e-^^t [I _ K) . (287) 

Zj 



If the impedance of the load is Z/,, then: 

Zl = j = Z^^. (288) 

Solving this equation for K (which gives the relative amplitude and phase of the "reflected" wave), we 
find: 

Zl/Z - 1 Zl- Z 

""-zU^i-z^- ^^''^ 

Note that if Zl = Z, there is no reflected wave: the termination is correctly matched to the characteristic 
impedance of the transmission line. 

11.3 "Lossy" transmission lines 

So far, we have assumed that the conductors in the transmission line have zero resistance, and are sep- 
arated by a perfect insulator. Usually, though, the conductors will have finite conductivity; and the 
insulator will have some finite resistance. To understand the impact that this has, we need to modify our 
transmission line model to include: 



a resistance per unit length R in series with the inductance; 

a conductance per unit length G in parallel with the capacitance. 



The modified transmission line is illustrated in Fig. 25 



The equations for the current and voltage are then: 



dV dl 

— = -L—-RI, (290) 
ox ot 
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d]_ 

dx 



-C 



dV 

~dt 



GV. 



(291) 



We can find solutions to the equations (290 1 and (291 1 for the voltage and current in the lossy transmis- 
sion line by considering the case that we propagate a wave with a single, well-defined frequency uj. In 
that case, we can replace each time derivative by a factor —iu:. The equations become: 



where 



dV 

dx 

dx 
L = L 



iujLI - RI 



iioCV - GV 



-L 



-G 



dV 

'dt' 



R 



and 



G = G 



G 



(292) 



(293) 



(294) 



The new equations ( 292 1 and (293 1 for the lossy transmission line look exactly like the original equations 



( 271| ) and (272i for a lossless transmission line, but with the capacitance G and inductance L replaced 
by (complex) quantities C and L. The imaginary parts of C and L characterise the losses in the lossy 
transmission line. 

Mathematically, we can solve the equations for a lossy transmission line in exactly the same way 
as we did for the lossless line. In particular, we find for the phase velocity: 

1 1 



(295) 



and for the impedance: 



Z 




L + i 



R 



(296) 



Since the impedance ( |296| ) is now a complex number, there will be a phase difference (given by the 
complex phase of the impedance) between the current and voltage in the transmission line. Note that 
the phase velocity ( 295| ) depends explicitly on the frequency. That means that a lossy transmission line 
will exhibit dispersion: waves of different frequencies will travel at different speeds, and a the shape of 
a wave "pulse" composed of different frequencies will change as it travels along the transmission line. 
Dispersion is one reason why it is important to keep losses in a transmission line as small as possible (for 
example, by using high-quality materials). The other reason is that in a lossy transmission line, the wave 
amplitude will attenuate, much like an electromagnetic wave propagating in a conductor. 

Recall that we can write the phase velocity: 



where k is the wave number appearing in the solution to the wave equation: 



and similarly for the current I. Using Eq. ( 295 1 for the phase velocity, we have 



k 




l + i 



R 
ujL 



l+i 
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(297) 



(298) 



(299) 



Let us assume R <^ loL (i.e. good conductivity along the transmission line) and G <C uiG (i.e. poor 
conductivity between the Unes); then we can make a Taylor series expansion, to find: 



k w ujVLG 



1 + 



2uj 



R G 



(300) 
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Fig. 26: Coaxial cable transmission line. 



Finally, we write: 



k = a + 1/3, 



and equate real and imaginary parts in equation ( |300| ) to give: 

Q ~ ujV LC, 

and: 



/3 



1 / R 



2 \Zi 







+ GZo], 



(301) 



(302) 



(303) 



where Zq = \/L/C is the impedance with R = G = (not to be confused with the impedance of free 
space). Note that since: 

V = Foe'e^^-'^*) = yoe"^'^e^("^-'^*\ (304) 
the value of a gives the wavelength A = 1/a, and the value of /? gives the attenuation length 6 = 1/(3. 



11 A Example: a coaxial cable 

A lossless transmission line has two key properties: the phase velocity v, and the characteristic impedance 
Z. These are given in terms of the inductance per unit length L, and the capacitance per unit length C : 



1 



z 



The problem, when designing or analysing a transmission line, is to calculate the values of L and C. 
These are determined by the geometry of the transmission line, and are calculated by solving Maxwell's 
equations. 

As an example, we consider a coaxial cable transmission line, consisting of a central wire of 
radius a, surrounded by a conducting "sheath" of internal radius d - see Fig.|26] The central wire and 
surrounding sheath are separated by a dielectric of permittivity e and permeability fi. Suppose that the 
central wire carries charge per unit length +A, and the surrounding sheath carries charge per unit length 
—A (so that the sheath is at zero potential). We can apply Maxwell's equation ([T]) with Gauss' theorem, 
to find that the electric field in the dielectric is given by: 



A 

27rer ' 



(305) 
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Fig. 27: Inductance in a coaxial cable. A change in the current / flowing in the cable will lead to a change in the 
magnetic flux through the shaded area; by Faraday's law, the change in flux will induce an electromotive force, 
which results in a difference between the voltages V2 and Vi. 



where r is the radial distance from the axis. The potential between the conductors is given by: 

V=l E-df=— In (-\ (306) 
Ja 27re \aj 

Hence the capacitance per unit length of the coaxial cable is: 

To find the inductance per unit length, we consider a length / of the cable. If the central wire 
carries a current /, then the magnetic field at a radius r from the axis is given by: 

\B\ = ^. (308) 



The flux through the shaded area shown in Fig. 27 is given by: 

^ = 1 j \B\dr = ^ In . (309) 

A change in the flux through the shaded area will (by Faraday's law) induce an electromotive force 
around the boundary of this area; assuming that the outer sheath of the coaxial cable is earthed, the 
change in voltage between two points in the cable separated by distance / will be: 

/\V = V2-Vi = ~. (310) 
at 

The change in the voltage can also be expressed in terms of the inductance per unit length of the cable: 

Ay = -/L^, (311) 
where / is the current flowing in the cable. Hence the inductance per unit length is given by: 



With the expression in Eq. (307 1 for the capacitance per unit length: 

C 



In {d/a) ' 

the phase velocity of waves along the coaxial cable is given by: 

1 1 



y/LC y/Jj^ 
and the characteristic impedance of the cable is given by: 



(313) 
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